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Abstract 

We consider the effects of quark masses to the perturbative thrust in e + e~ annihila- 
tion. In particular we show that perturbative power corrections resulting from non-zero 
quark masses considerably alters the size of the non-perturbative power corrections and 
consequently, significantly changes the fitted value of a s . 
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One of the cleanest signatures of perturbative QCD comes from jet cross sections in e + e _ 
annihilation. In such processes, it is possible to define infra-red safe event shape variables 
which can be calculated order by order in perturbative QCD and compared subsequently with 
experiment. However in order to carry out these comparisons, a method has to be evolved to 
parametrise non-perturbative effects which though expected to be small at present Q 2 values 
at LEP, actually turn out to be substantial (~ 25%) even at Q ~ mz- One of the reasons for 
this is that these non-perturbative effects are actually suppressed by a single power of Q rather 
than Q 2 . In addition, it is also possible that these power corrections could be comparable to 
0(a 2 ) at present LEP energies. 

In order to address these issues the Milan group of Dokshitzer et al. 1J drawing on the 
earlier work of Webber [2 , Korchemsky and Sterman [3| and others, presented a systematic 
approach for handling power corrections using perturbation theory. Very briefly, they studied 
the consequences of assuming that a s has a low energy effective form which does not grow at 
low scales but has an infra-red regular form. The moments of a s are integrated only over the 
infra-red region. Various non-perturbative parameters are then parametrised and the form and 
magnitude of power corrections are determined. 

However before one uses the approach of the Milan group in order to get a handle on power 
corrections and subsequently determine a s by a fit to the data, it is important to isolate power 
corrections coming from a purely perturbative region. The Milan approach neglects the masses 
of all the quarks but instead uses a gluon mass as a 'trigger' to differentiate the perturbative 
from the non-perturbative region. We find however that the masses of the quarks, particularly 
the c and the b quarks, even at present LEP energies, can contribute significantly (of the order 
of about 25%). In fact, if we go beyond the top quark threshold (which is expected, perhaps in 
the future NLC) the perturbative contribution to power contributions due the top quark mass 
is even larger. We will have more to say on this later in the paper. 

In this paper, we consider the example of one such event shape variable - the thrust - and 
show the significance of the effect of quark masses which need to be folded in before estimating 
the non-perturbative contribution to power corrections. We present explicit expressions to 
0(a s ) of quark mass corrections expanded to 0(m). We also show the effect of keeping the 
full mass contribution to 0(a s ) which unfortunately does not have a simple analytic form like 
the former and needs to be calculated numerically. Using these expressions we then fold in the 
power corrections of the Milan type and use this full expression to estimate both a and a s 
and compare it with estimates that exist in the literature without taking quark masses into 
account. 

The first paper which calculated the effect of quark masses to 0(a s ) was published about 16 
years ago by one of the authors 0]. For completeness, in what follows, we quote those results 
from that paper which we need for our analysis here. The thrust, as defined traditionally, is 
given by 

_ n maxJ2 ith ( Pi ■ n) 
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where the denominator runs over all observed particles and the numerator runs over all particles 
in a hemisphere, n is a unit vector chosen in a direction that maximises the numerator and 
defines the jet axis. 

While this definition is appropriate for all massless particles, to include mass effects in the 
definition of the thrust, we modify the above definition slightly and write 
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where W 2 = s. Of course the denominator equals J2i \Pi\ when all the particles are massless. 
This normalisation with the total energy is also what is used by the Milan group in their analysis 
though in their case the massive gluon eventually decays into massless quarks and gluons. 
For a three particle final state, the thrust, as we define it, is given by 



T = max 
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where X{ = 2Ei/W, being the energy of the zth particle in the final state in the cm. frame and 
£ = 4m 2 /W 2 , m being the mass of the quarks. Note that in the two-jet limit T = Tq = y/1 — £. 
The average value of the thrust is defined by 



<T >- 
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The numerator of the above is given up to 0(a s ) and to 0{£) by (cr = [Ama 2 / s)e\ is the total 
cross section for e + e~ — > q^i) [I] 
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where Li 2 (x) is the dilogarithm function. In the £ — > limit this gives, for the average thrust, 



<T>=1 + -Ji 
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which works out to, for the perturbative thrust in the massless limit, 



< 1 - T >= 1.05—, 
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as quoted in numerous places in the literature. 

Several points here are worthy of note. The leading term in the 0(m) expansion above 
is £ x / 2 . Thus the leading mass correction goes as 1/Q. To the best of our knowledge, this 
fact was noticed for the first time in and subsequently in [2] and |1J who have traced it to 
appear from the soft phase space boundary. We would like to stress that this 1/Q behaviour 
is a pure perturbative higher twist effect to the thrust and not related to any non-perturbative 
contribution. Thus, it seems clear, that the coefficient of 1/Q in the full expression for the 
thrust would include contributions both from the perturbative as well as the non-perturbative 
sectors. This aspect will become more quantitative, when we do our fits later. 

The second point to note is a calculational one. Since, in the two jet limit, the thrust is 
equal to To = y/1 — £, in order to make the virtual contributions vanish we need to calculate, 
not as in the usual case < 1 — T >, but < To — T >. It is then a trivial matter to add a term 
< 1 — Tq > to obtain < 1 — T > to compare with experiment. 
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Figure 1: Effect of mass correction and higher order in the mean thrust value. The solid and 
the dashed lines are respectively the first and second order calculations of < 1 — T > without 
mass corrections. The dashed-dotted and dotted lines are complete and approximate (0(m)) 
first order calculations with quark mass effects. 



In order to compare with experiment, however, and to redo the fits for a s and «o we have 
used not only the 0(m) contribution above but also the full massive contribution to 0(a s ) albeit 
evaluated numerically. In addition we have also compared the 0(a^) massless corrections to 
the 0(a s ) massive correction to try and estimate how much of the 1/Q corrections can be 
mimicked by higher orders in the coupling constant. 

Figure ^ shows < 1 — T > as a function of centre of mass energy computed with a s (mz) = 
0.12. As one sees from the curve the contribution of the second order terms are large over the 
entire energy region (~ 55% at Q = 12 GeV going down to 33% at 200 GeV). On the other 
hand the effect of quark masses, evaluated only to first order in a s , is even larger at low centre 
of mass energy (~ 76% at Q = 12 GeV). This is clearly a 1/Q power law effect and hence dies 
off faster, becoming 2.5% at Q ~ m z and negligible at 200 GeV. It is clear from the figure that 
one needs to take the full massive correction rather than the 0(m) contribution, because it 
accounts only 60% (30%) of the mass correction at 20 GeV (12 GeV). 

In order to compare the theoretical predictions with the measurements done at different cen- 
tre of mass energies [Sj at PETRA, PEP, TRISTAN, SLC and LEP, we add the non-perturbative 
contribution a la the Milan group [T to the perturbative contribution. In this paper, we use 
only the 0(a s ) calculation of < 1 — T > and a more detailed comparison with a 0(a^) calcu- 
lation is under preparation [6 . We will have more to say on this later. The non-perturbative 
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contribution, as is well known, is given by an additive contribution < 1 — T > pow : 

« b 2 (Q) 
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where \ij is an infra-red matching scale (taken as 2 GeV), K = (67/18 — iy 2 /6)-Ca — 5N f /9 
and M. is the Milan factor (determined to be 1.49) [Zj. 
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Figure 2: Results of the fit of < 1 — T > to the first order calculations with and without quark 
mass effects. The points indicate measured < 1 — T > values from different experiments and 
the solid and dashed lines are fits without and with mass corrections. 

Figure |2] shows the experimental values of < 1 — T > together with the two fits which use 
respectively the massless and the massive forms (both to 0(a s ) for the perturbative contribu- 
tion). These fits have been carried out with two free parameters a s (mz) and «o- Both massless 
and massive formulation of the perturbative component give reasonable fits to the data with 
X 2 of 90.2 and 69.2 respectively for 48 degrees of freedom corresponding to confidence levels of 
0.22 x 10 -3 and 0.24 x 10 _1 . However, they do differ in the final values of a s (mz) and a as 
can be seen in Tabled In these fits the scale parameter is chosen to be 1.0. 

The errors quoted in the Table Q are experimental errors obtained from the minimisation 
procedure. We can also estimate the theoretical uncertainties on these quantities by varying the 
scale parameter. If we vary the scale parameter between 0.5 and 2.0, we obtain uncertainties 
in a s and a Q to be ±0.010 and ±0.12 respectively. The value of a s (m z ), obtained from the fits, 
when quark mass effects are included or ignored, differ by 0.008 which is much larger than the 
experimental uncertainty of about .001 on the a s value and comparable in fact to the theoretical 
uncertainty. 
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Fit Type 


a 


a s (m z ) 


Massless quarks 
Massive quarks 


0.7931 ±0.0066 
0.7385 ±0.0065 


0.1533 ±0.0015 
0.1612 ±0.0014 



Table 1: Results of the fits of the < 1 — T > distribution to perturbative and power law terms 
when the quark mass effects in the perturbative term is ignored or included 

It is thus clear from the preceding analysis that an estimate of the power corrections due to 
the non-zero masses of the quarks is crucial in getting better and more realistic estimates on 
the strong coupling constant and indeed, in general, on power corrections. The next obvious 
step would be to calculate mass corrections to 0(a^). Some results in this direction have been 
obtained by Nason and Oleari [Hj which could be used to carry out a similar analysis to the one 
presented above. We are, at present, in the process of extending our analysis to second order 
in the strong coupling using the results of jH]. 

In various projected Linear Collider scenarios (like, for example, the NLC) energies upwards 
of 500 GeV are expected. In such a region, the effect of the top quark would be dramatic and 
significant. The combination of the large mass of the top quark and a charge squared of 4/9 
implies that the usual massless expressions for the thrust would not work. We have estimated 
that the difference between choosing a massless formula for describing the thrust beyond the 
top quark threshold and using the (more appropriate) massive formula changes the value of 
the thrust by about a factor of 5 near the threshold. Most of this contribution comes, in fact, 
from the top quark mass. In the table below we give an estimate of the change that would 
occur between choosing all quarks massless and massive above the top quark threshold. It is 
obvious that the effect is spectacularly large, particularly near the threshold. Mass effects in the 
resummation of event shape variables are also expected to be significant and this is presently 
being studied. 



Q (GeV) 


a s 


< 1 — T > (massless) 


< 1 — T > (Massive) 


360 


.0995 


.0323 


.1605 


500 


.0957 


.0311 


.0994 


1000 


.0888 


.0289 


.0435 



Table 2: Difference between choosing massless and massive quarks above the top quark thresh- 
old 

Thus, it is imperative that in order that reliable estimates be made of the thrust at these 
energies, we have available, calculations to higher orders in a s of e + e~ scattering with massive 
quarks in the final state. This would also give us a handle on the relative magnitudes of power 
corrections to the thrust to a particular order in a s and the magnitude of the next order term 
in a s [HE!- For example NNLO effects might be capable of mimicking the 1/Q behavior. Mass 
effects in the resummation of event shape variables are also expected to be significant and this 
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is presently being studied. 

Acknowledgements: One of us (R.B.) would like to thank George Sterman for many illumi- 
nating discussions. This project was started at the Sixth Workshop on High Energy Physics 
Phenomenology (WHEPP-6), held at Chennai, India in January 2000 and we would like to 
thank all the funding agencies which made this workshop possible. 



References 

[1] Yu. L. Dokshitzer and B. R. Webber, Phys. Lett. B352 (1995) 451; 

Yu. L. Dokshitzer, G. Marchesini and B. R. Webber, Nucl. Phys. B469 (1996) 93. 

[2] B. R. Webber, Phys. Lett. B339 (1994) 148. 

[3] G. P. Korchemsky and G. Sterman, Nucl. Phys. B437 (1995) 415. 
[4] R. Basu, Phys. Rev. D29 (1984) 2642. 

[5] ALEPH Collaboration, D. Bukusulic et al, Z. Physik C55 (1992) 209. 

ALEPH Collaboration, D. Bukusulic et al, ALEPH 98-025 (1998), contributed to 199£ 
winter conferences. 

AMY Collaboration, Y. K. Li et al, Phys. Rev. D41 (1990) 2675. 
CELLO Collaboration, H. J. Behrend et al, Z. Physik C44 (1989) 63. 
DELPHI Collaboration, P. Abreu et al., Phys. Lett. B456 (1999) 322. 
JADE Collaboration, A. Morivilla Fernandez et al., Eur. Phys. J. CI (1998) 1. 
L3 Collaboration, M. Acciari et al., CERN-EP/2000-064 (2000). 
MARK-II Collaboration, A. Peterson et al, Phys. Rev. D37 (1988) 1. 
MARK-II Collaboration, A. Abrams et al, Phys. Rev. Lett. 63 (1989) 1558. 
MARK-J Collaboration, D. P. Barber et al, Phys. Rev. Lett. 43 (1979) 902. 
OPAL Collaboration, G. Abbiendi et al, CERN-EP/1999-178 (1999). 
PLUTO Collaboration, Ch. Berger et al, Z. Physik C12 (1982) 297. 
TASSO Collaboration, W. Braunschweig et al, Z. Phys. C47 (1990) 187. 

[6] Under preparation. 

[7] Yu. L. Dokshitzer, |hep-ph/9911299| 



[8] P. Nason and C. Oleari, Nucl. Phys. B521 (1998) 237. 

[9] G. Sterman, talk at the QCD-Euronet workshop, Florence, Italy, October 1999. 



6 



